Dynamic phase transition in the conversion of B-DNA to Z-DNA 
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The long time dynamics of the conformational transition from B-DNA to Z-DNA is shown to 
undergo a dynamic phase transition. We obtained the dynamic phase diagram for the stability of 
the front separating B and Z. The instability in this front results in two split fronts moving with 
different velocities. Hence, depending on the system parameters a denatured state may develop 
dynamically eventhough it is thermodynamically forbidden. This resolves the current controversies 
on the transition mechanism of the B-DNA to Z-DNA. 
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I. INTRODUCTION 

The most common form of DNA found under normal 
physiological low salt conditions is the B-DNA, the well 
known double helix with a right handed helicity. Quite 
surprisingly, the first DNA structure to be solved by X- 
ray crystallography turned out to be a left handed zig- 
zag form called the Z-DNA [^Q. This Z-DNA can be 
stabilized in vitro in presence of high salt concentration, 
cations or negative super-coiling. Although the Z form is 
transient in vivo due to the lack of friendly environment, 
still the B-Z conversion is relevant in poxviruses and 
Alzeimer's disease^. Apart from the inversion of helic- 
ity, the Z-DNA has a repeat unit of 2 base pairs compared 
to one for B. Thus a conformational transformation from 
the B to the Z form takes place as, e.g., the ionic concen- 
tration or super-coiling is changed. The B-Z transition 
is first order in nature [H, IB-El- 

As the base pairs and a subset of backbone sugar 
rings have to flip to execute the B-Z transition involving 
changes in helical chirality, the dynamics offers intriguing 
possibilities Only recently methods have been devel- 
oped to explore the dynamics in single DNA as opposed 
to earlier studies in solutions, though with conflicting re- 
sults. In Ref. [l(J the B-Z conformational transformation 
for a short 15 base pair GT (non- Watson-Crick wobble 
base pair) DNA wrapped on a single walled carbon nan- 
otubc was monitored as a function of time by the addition 
of counter-ions. The nanotube helped in identifying the 
phases via accurate measurements of the band-gap in a 
simpler geometry. This transition is completely reversible 
and is thermodynamically identical to the transition seen 
in the absence of the nanotube. The results seem to 
indicate the formation of a denatured DNA during the 
transformation, eventhough a denatured state under the 
experimental conditions is not possible thermodynami- 
cally. A different single molecule experiment studied the 
transition under a tension and negative super-helicity by 
combining FRET with magnetic tweezers [llj. This ex- 
periment on an effectively (GC)n DNA (i.e. 22 bases) 
seems to favour a single interface between B and Z with- 
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out any denatured bubble. 

One can characterize the B-Z transformation by a 
growth of a suitable domain over the bulk of DNA. In any 
such scenario, the B-Z interface, the separator between 
the two chiral phases, plays an important role. The equi- 
librium interface has been characterized structurally and 
from other studies. The structure of a short oligomer 
in presence of a Z-DNA binding protein at 2.6A reso- 
lution indicates broken base pairs separating the B and 
the Z phases. The protein acting as an external source 
inducing the transition is expected to produce a sharp 
interface fl2l|. A very ingenious way of studying the in- 
terface is to use mirror DNAfpjlldj]. though it cannot be 
used for chirality changing transition. Intcrfacial stud- 
ies and melting of short B-B* oligomers, where B* is 
the enantiomer of B, show that the junction mimics the 
B-Z junction, and the interface broadens as the melting 
point is reached. In contrast to these equilibrium cases, 
the nature of the interface during the transition depends 
on the nature of the transition mechanism0, H|. Sev- 
eral such schemes are in vogue and discussed in detail 
in Ref. @- The two main competing hypothesis for the 
B to Z transition mechanism are the following. It is ei- 
ther via (1) the base pair separation followed by base 
pair flipping [l[, or (2) the base pair flipping without any 
base pair separation^. In the first case there could be 
a denatured intermediate state while in the second there 
could a Z type but following the standard Watson-Crick 
base pairing (ZWC-DNA). 

The apparently contradictory results from the two sin- 
gle molecule experiments and the controversy associated 
with the B-Z transition mechanism in general, motivate 
us to study a coarse-grained thermodynamic model for 
the dynamics. This implicitly requires infinitely long 
chains, since small length DNAs or oligo-nucleotides do 
not show any proper thermodynamic phase transition. 
Inducing the B to Z transition is tantamount to a lower- 
ing of the free energy of Z compared to B making Z the 
most preferable state. We in our formulation use the sim- 
plification of the single molecule experiments to restrict 
the geometry to one-dimension only. The boundaries of 
the long chain are maintained in the two states so that the 
new structure develops from one side. In such a problem 
the dynamics of the transition produces a steady state 
with uniformly moving front (or fronts). An investiga- 
tion of the various types of fronts would help us answer 
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the question of any dynamic generation of thermodynam- 
ically forbidden state. 
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FIG. 1: (a)The Landau function F(<j>) as a function of (f>. In all 
cases, <j> = 1 is the stable state, Z-DNA, (f> — — 1 represents an 
unstable (dashed line) or metastable (solid and dotted lines) 
state, B-DNA while (f> = is a quadratically unstable (solid 
and dashed lines) or metastable (dotted line) state, denatured 
state. The three cases I, II and III in the text correspond to 
dotted, solid and dashed lines. (b)Potential V(U) = —F(U) 
for the particle-on-a-hill analogy. 



II. MODEL 

Our model consists of three states B, the denatured 
state and Z to be represented by a parameter <p = 
— 1,0,1. The space time coordinates z,t are taken to 
be continuous. It is a one dimensional problem where 
4>(z, t) describes the state of the coarse-grained base-pair 
at index z along the DNA. For the B-Z transition, we take 
6 = — 1(B state) to be unstable (or metastable) which is 
getting invaded by the stable state at cf> = 1(Z state). 
We study this phenomenon through a Landau free en- 
ergy F(<p) taken as a sixth order polynomial with the 
coefficients chosen to have extrema at <fi = 0, ±1. This is 
ensured by choosing the thermodynamic force f(4>) as 



m = 



dF{cp) 



+ a)(4> -/?)(! -«£)(! + 0), (1) 



where a, /3 > are constants, whose values are system 
specific. Needless to say, the relative stability of the three 
phases can be adjusted by a, ft. The Landau Ginzburg 
free energy is taken as 
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(2) 



where D > is the elastic constant. D-term allows inho- 
mogeneity, e.g., at the interface between two phases. The 
three homogeneous phases are given by the minima of the 
Landau free energy F(<f>). The dynamics is governed by 
the non linear diffusion equation 



d(f) d 2 (j) 



(3) 



derived from Eq. [5] in the overdamped limit. The friction 
coefficient has been absorbed in the definition of time t. 



The geometry to be considered is such that the B state 
is on one side and the Z state on the other with the front 
moving towards the unstable state. For the B-Z case, 
this is ensured by the boundary conditions 
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for Eq. [3] for all time. A few other boundary conditions 
are considered too. The three generic cases obtained by 
fixing a and (3 are the following(see Fig. [lji) 

• Case I : While quenching to the stable state, Z, 
state B remains in a metastable state while the de- 
natured state cf> — is also metastable. Since the 
barriers are somewhere in between (f> = — 1 and 
4> = 1, we have < a, j3 < 1. 

• Case II : The metastable state(B-DNA) sees a bar- 
rier somewhere inbetween —1 to 0, while the de- 
natured state is quadratically unstable state. This 
case is for < a < 1, and (3 = 0. 

• Case III : Unstable B state quenched into stable 
Z while the denatured state remains in a quadrat- 
ically unstable state (i.e., without facing any bar- 
rier). This happens when a > 1 and /3 = 0. 

To be noted that cases I and II are similar to the free 
energy landscape obtained in Ref. [ll| as the potential of 
the mean force obtained from molecular dynamics. 

The diffusive like term in Eq. [3] coming from the elastic 
part of Eq. [2] tends to smoothen out any inhomogeneity 
while the driving force f(4>) tends to favour the stable 
state whenever there is any inhomogeneity. The com- 
bined effect of the diffusion like spreading and the selec- 
tion of one phase by the drive leads to a steady state 
where the interface shows a uniform motion and takes a 
shape which is not necessarily the equilibrium shape fl6|. 
Based on the Fishcr-Kolmogorov(F-K) idea, the traveling 
wave solution 4>{z,t) = U{z — vt) can be used to rewrite 
Eq.GHas 



dU *,TT, n 



(r = z — vt), 



(4) 



where v the velocity of the front is to be determined. The 
interface which we are studying is between <p = +1 and 
<j) = — 1 states. Eq. [l]can be interpreted as the motion 
of a particle moving in a potential V = —F(U)(Fig. [Tb) 
starting at the hill at U = +1 at time r = — oo just 
reaching the other hill at U = — 1 at time t = +oo los- 
ing energy due to "friction" v. For a given potential, 
such a motion is possible only for particular values of 
v and that velocity is the selected velocity of the front. 
However, it is also possible that the particle spends an 
infinite amount of time in the intermediate state so that 
the descent from U = +1 to U = and the descent 
from U = to U = —1 are independent requiring two 
different friction coefficients. The physical picture that 
emerges is that the stable state moves towards the un- 
stable state, and the propagating front will have a time 
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the limit of stability of the BZ front(from below). 



FIG. 2: (a)Plot of velocity vs a for a fixed value of p = 
0.45. Three velocities meet at a common point at a c (j3). The 
remaining three figures(0 vs z) represent the time evolution 
of the front(or fronts). (b)A single front for a = 0.6 < a c (j3). 
(c) A single front for a = 0.7 near a c (/?) with a signature of the 
width widening but no "0" phase. (d)For a — 0.72 > a c 
single front splits into two fronts. 



independent shape and a constant velocity v. However 
in some situations, the initial big front separating the 
two phases = ±1 splits into two, one front between 
= — 1 and = 0, while the other one between = 
and = 1. The two smaller fronts move with different 
shapes and speeds w_io,woi- The = state may then 
get dynamically generated. Consequently one may see 
the development of the denatured state. The less prefer- 
able state will eventually be devoured by the stable state 
completing the transition from B to Z-DNA. 



III. DYNAMIC PHASE 
DIAGRAM-NUMERICAL AND PERTURB ATIVE 
APPROACH 

The velocity of the front has been determined by nu- 
merical analysis for different boundary conditions like 

(a) (j)(—oo,t) — 1, 0(oo,t) = —1 for the B-Z front, 

(b) 0(— oo,t) = 1, 0(oo,i) = for a front between Z and 
the denatured state, (c) </>(— oo,t) = 0, 0(oo,i) = —1 
for a front between B and the denatured state. The 
initial(t = 0) interface of width w is located at z = zq and 
a Crank Nicolson method is used to evolve the nonlinear 
diffusion equation. Once a steady state is reached, the 
velocity is determined by locating the positions at which 
4> = ±.5, and = as appropriate. In the case of the 
split front, only the velocity vq\ can be obtained by the 
F-K analysis but not in general. 

The dependence of the velocities on a for a fixed (3 is 
shown in Fig. [5^,. We see that three fronts move with 
different velocities for a < a c ((3) with «oi > u_n > 
v_io- All these velocities are same at a = a c (/3). For 
a > a c (/3), the BZ front splits into two fronts and the 
denatured state grows with time as (v_io — voi)t. It is 
straightforward to see that no stable front between ±1 
can exist if V-iq < vqi- Also the V-n curve ends at 
u c {f3) and has no continuation for a > a c {(3). This in- 
dicates that a c (/3) is a singular point. The numerically 
determined a c (J3) vs (3 line is shown in Fig.[3J This is the 
phase diagram for dynamics with the phase boundary as 




FIG. 3: Dynamic phases in a plot of a vs /3, the bound- 
ary (solid line) being given by a = a c (/3). In the region below 
the boundary line, a single front between —1 to l(big front) 
propagates without splitting. In the region above the bound- 
ary line the front between —1 to +1 splits into two (small) 
fronts. Z, "0" and B are the stable states in regions 1, 2 and 
3 respectively. The dotted line corresponds to Vqi = 0, while 
the dash-dotted line to f-io = 0. The split fronts move away 
from each other in region 2, both towards right in I and both 
towards left in 3, as per the chosen boundary conditions. The 
big front has zero velocity on the a — j3 line and the diagram 
is symmetric around this line. Point q represents the equilib- 
rium point, where three states have the same free energy. 

The phase diagram can be confirmed by considering a 
few special cases. For a = (3, the free energies of B and 
Z are same and the BZ front should have zero velocity. 
The point a = f3 = corresponds to the equilibrium 
situation, for which all the three fronts are static, and 
therefore the condition to be on the phase boundary is 
trivially satisfied. This point is denoted by q in Fig. |3] 
Along the a = /? line for a < a c (/3), v i,V- W ^ with 
state +1 or —1 invading 0. In contrast in region 2, along 
the same a = (3 line, "0" is the stable state and it invades 
both ±1 states. In region 2 above the dotted line, ob- 
tained by equating F(l) = F(0) (Eq. {TJ), the "0" state 
grows with the two fronts moving away from each other, 
but below that dotted line in region 1 the Z state grows 
though the fronts move in the same direction(towards 
right). The Z B symmetry in our choice of the free 
energy mandates a symmetric phase diagram across the 
a = j3 line with the fronts moving towards left in region 
3. 

For a,j3 close to a = /3 = -A=, a pcrturbative 



analysis [17] can be done to determine the velocity, which 
is now a small parameter. By writing, to first order in v, 



U(z- vt) w U 0) + V t Uq (z) 



(5) 



where prime denotes a derivative, v can be determined 
to first order in free energy difference if Uq is known. 
In the equilibrium situation, there is a Goldstone like 
zero-energy mode, because, the interface can be placed 
anywhere or shifted along z without any cost of energy. 
We therefore take Uq{z) as centered around an arbitrarily 
chosen origin. The static solution satisfies, 



F([/ ) = [/ 2 ([/ 2 -l)' 



(6) 
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With a first order correction, the velocities are 



v 



<j 



i j 



(7) 



where i,j = 0, ±1, and the free energy differences are 
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(8) 
(9) 
(10) 



At this perturbative regime, by equating the velocities, 
we find that around a = /? = the slope of the crit- 
ical line is —1, which is consistent with the numerically 
determined boundary shown in Fig. [31 Moreover we also 
find the phase boundary to deviate very slightly from a 
straight line over the range shown there. There is a devi- 
ation from linearity beyond that but the numerical error 
becomes larger. We next study the behavior of the width 
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FIG. 4: (a) Time evolutions of the width (in arbitrary units) 
of the front are shown for different values of a keeping /3 = 
0.4 fixed. The initial front had a stretch of "0" phase which 
decays for a < a c (fi) but grows linearly for a > a c (/?). The 
timescale to reach saturation increases as a — > a c (/?) — . For 
this case a c (/3) sa 0.752. (b) Log-Log plot of width versus 
time scale for a < a c (/3). Two solid line slopes are shown far 
from and near a c (/3) . 

of the interface and of the appropriate timescale for the 
dynamics. For the special case of /3 = 0.5 as a — > a r (/3 ) 
the divergence of the width has been noted in Ref. [181 ] . 
At a = P = y-, any length of "0" domain can be inserted 

and therefore the width of the BZ interface at the limit of 
stability is infinity. On the split-front side (Fig. [5Ji), the 
width increases linearly with time as W = (vqi — V—i t o)t 
(Fig. 0Ji for a = 0.75205). While, on the other side of 
the phase boundary the single front(Fig. [5}}) has a finite 
width(Fig. [4ji for a — 0.745). Close to the phase bound- 
ary though a deformation of the moving front is visible 
around <fr = 0(Fig.[3b), but width saturates at large time 
(Fig. 0k for a = 0.75175) without any appearance of the 
denatured phase. Hence scaling forms are expected as 



Fig. |4^ shows the time evolution of the width of an inter- 
face for various a at a fixed (3, where the instantaneous 
width W of the interface at time t is obtained as 



W 2 = < z 2 > - < z > 2 , where 

2 



< z" > = 



dz 



d<t>(z,t) 
dz 



(11) 

(12) 



Another way to characterize the width would be to look 



at the slope of the profile 



d<t>(z) 



dz 



which is re- 



W ~| a-a c ((/3) |" 



W z . 



lated to the inverse of W and also shows the scaling with 
characteristic dynamic exponent. We started with an in- 
terface that has an insertion of the "0" state and the 
width monitors the decay or the growth of the "0" state. 
The width saturates exponentially for a < a c albeit 
slowly near a — > a c — , while a linear growth is observed 
for a > a c (j3). Time here refers to the discretized time in 
the Crank-Nicolson approach. By fitting an exponential 
to the time evolution of W, the characteristic time scale 
was determined, for a < a c (/?) . The exponent (i is found 
to be rather small, not inconsistent with the logarithmic 
growth observed in Ref. [HI]. Fig. HJd shows the log-log 
plot of r vs W indicating a value of z within 3.0 to 4.0. 
However for better accuracy one requires a large system 
and long time observation as well. The divergences of 
W and r with scaling establish the critical nature of the 
a = a c (/3) line. 

Despite the immense success in probing the various 
phases of DNA by single molecular manipulation tech- 
niques, interfaces have not been explored thoroughly. We 
hope our results will motivate direct studies of interfaces 
of DNA, especially their stability. Even on the theoretical 
front, it remains to be seen if all atom molecular dynam- 
ics simulations that have been successful [HI, [JJ| [2(| in 
seeing various phases, can be used to monitor the dynam- 
ics of interfaces, B-Z in particular, under given boundary 
conditions. 



IV. CONCLUSION 

The conformational transition from B-DNA to Z-DNA 
has been studied via wave-front propagation. The dy- 
namic phase diagram for the steady state is obtained in 
the a-/3 plane, where a, (3 characterize the relative sta- 
bility of the phases, by the critical value a c for differ- 
ent values of /3. The phase boundary in the a-/3 plane 
has been determined and corroborated by a perturbation 
analysis. The dynamic transition is associated with di- 
verging length and time scales and has its own dynamic 
exponent. On one side of the phase boundary the dy- 
namics involves propagation of one B-Z interface with a 
uniform speed, while on the other phase such an interface 
is unstable leading to the formation of the thermodynam- 
ically forbidden denatured state. This in turn, suggests 
that there is no unique mechanism for the B-Z dynamics 
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and it is possible to switch from one type to other by in experiments is that the two cases, namely nanotubc 
tuning the parameters. A resolution of the controversy and magnetic tweezers are on the two sides of the phase 

boundary. 
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